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Abstract:  This  paper  studies  robustness  of  multivariable  systems  with  parametric  uncer¬ 
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1  Introduction 

Motivated  by  the  seminal  theorem  of  Kharitonov  on  robust  stability  of  interval  polynomials[l, 
2],  a  number  of  papers  on  robustness  analysis  of  uncertain  systems  have  been  published  in  the  past  few 
years[3,  4,  5,  6,  7,  8,  9,  10].  Kharitonov’s  theorem  states  that  the  Hurwitz  stability  of  the  real  (or  complex) 
interval  polynomial  family  can  be  guaranteed  by  the  Hurwitz  stability  of  four  (or  eight)  prescribed  critical 
vertex  polynomials  in  this  family.  This  result  is  significant  since  it  reduces  checking  stability  of  infinitely 
many  polynomials  to  checking  stability  of  finitely  many  polynomials,  and  the  number  of  critical  vertex 
polynomials  need  to  be  checked  is  independent  of  the  order  of  the  polynomial  family.  An  important 
extension  of  Kharitonov’s  theorem  is  the  edge  theorem  discovered  by  Bartlett,  Hollot  and  Huang[4].  The 
edge  theorem  states  that  the  stability  of  a  polytope  of  polynomials  can  be  guaranteed  by  the  stability 
of  its  one-dimensional  exposed  edge  polynomials.  The  significance  of  the  edge  theorem  is  that  it  allows 
some  (affine)  dependency  among  polynomial  coefficients,  and  applies  to  more  general  stability  regions, 
e.g.,  unit  circle,  left  sector,  shifted  half  plane,  hyperbola  region,  etc.  When  the  dependency  among 
polynomial  coefficients  is  nonlinear,  however,  Ackermann  shows  that  checking  a  subset  of  a  polynomial 
family  generally  can  not  guarantee  the  stability  of  the  entire  family[ll,  12]. 

For  Hurwitz  stability  of  interval  matrices,  Bialas  ’proved’  that  in  order  to  guarantee  robust  stability, 
it  suffices  to  check  all  vertex  matrices [13].  Later,  it  was  shown  by  Barmish  that  Bialas’  result  was 
incorrect [14].  Kokame  and  Mori  eastblished  a  Kharitonov- like  result  on  robust  Hurwitz  stability  of 
interval  polynomial  matrices [15],  and  Kamal  and  Dahleh  established  some  robust  stability  criteria  for 
MIMO  systems  with  fixed  controllers  and  uncertain  plants  [16]. 

In  this  paper,  we  will  study  robustness  of  a  class  of  MIMO  systems  with  their  transfer  function 
matrices  described  by 

{/  an(s)  ...  ain(s)  \ 

I  •  •  •  .  I  •  aij(s )  €  Aj,j(s) 

\  Q“nl  (s)  Unn(s)  J 

Ay(s)=conv  {bli:j(s), brriij(s)} 
where  m  is  a  given  positive  integer. 

2  Preliminaries 

Definition  1  A  polynomial  matrix  is  a  matrix  with  all  its  entries  being  polynomials. 

Definition  2  Suppose  I?  is  a  simply-connected  region  in  the  complex  plane.  If  all  the  roots  of  the 
determinant  of  a  polynomial  matrix  lie  within  D,  then  this  polynomial  matrix  is  called  D-stable.  A  set 
of  polynomial  matrices  is  called  robustly  D-stable,  if  every  polynomial  matrix  in  this  set  is  D-stable. 
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Definition  3  Suppose  . . . ,  /m(s)  are  m  given  polynomials,  the  set 

m  m 

Y  A ifi(s)  :  A i  >0,  Y,  Aj  =  1 

i—1  i= 1 

is  called  the  polynomial  polytope  generated  by  fi(s), . . . ,  /m(s),  denoted  as  conv  {/i(s), . . . ,  /m(s)}. 

Definition  4  The  polynomial  /(s)  =  a0  +  aiS  +  . . .  +  ansn,  with  a,  €  [af ,  a\J]  is  called  an  interval 
polynomial. 

Definition  5  Sn  is  the  set  of  all  bisections  from  {1, . . . ,  n}  to  {1, . . . ,  n}. 

Definition  6  The  vertex  set  and  edge  set  of  Atj  (s)  are 

E%j  ( S ^  f  i  ?  ( )  ;  •  *  •  >  b  'f  f  l  ij  (  s  ])  | 

Eij(s)  =  {Xbrij(s)  +  (1  -  X)btij(s),  As  [0,1],  r,t  G  {1, . . .  ,m}} 

respectively. 


Definition  7 


>)  \  (aij(s))nxn  •  G'tj(s) 

<t£S„  ^ 


G  Eij(s)  if  i  =  a(j) 
G  Kij(s)  if  i  ^  a(j) 


(2) 


Lemma  1  (Edge  Theorem[4])  Suppose  T  C  C  is  a  simply-connected  region,  Q  is  a  polynomial 
polytope  without  degree  dropping.  Then,  f l  is  T-stable  if  and  only  if  all  the  edges  of  f l  are  T-stable. 

Lemma  2  Suppose  A(s)  is  a  given  n  x  (n  —  1)  polynomial  matrix.  Then 

on(s) 


i  Ms)  I  : 

Onl(s) 

<G>  for  alii  =  1, . . . ,  n, 


an(s)  G  An (s) , 
i  =  1, ...  ,7i 

an(s) 


is  robustly  D-stable 


A(s)  1  :  TiU  G  E,j(sj  \t\  \  is  robustly  stable. 


ani(s) 

Proof:  Necessity  is  obvious,  since  the  later  is  a  subset  of  the  former. 
Sufficiency:  For  any  a,ji(s)  G  An  (s),  the  corresponding  matrix  is 

an(s) 


T(s)  = 


Ms) 


«nl(s) 


By  Laplace  Formula,  we  can  expand  the  determinant  of  T(s)  along  its  first  column.  Then,  by  convexity 
and  by  Lemma  1,  we  know  that  T(s)  is  robustly  D-stable. 

Lemma  3  Suppose  B(s)  is  a  given  (n  —  1)  x  polynomial  matrix.  *  stands  for  fixed  entries  in  a 
matrix.  Then 


*  a±i(s)  *  aij(s)  * 
B(s) 

*  aij(s)  *  aij(s)  * 

B(s) 


is  robustly  D-stable 


Oi,(s)  G  xfii(s) 
aij(s)  G  A!j{s) 

ai,(s)  x  ai^s)  G 

(iFii(s)  x  Eij(s))  U  {Eu(s)  x  Kij(s)) 


is  robustly  D-stable. 


Proof:  the  proof  is  analogous  to  the  proof  of  Lemma  2,  except  that  the  Laplace  expansion  is  carried 
out  along  the  row  instead  of  the  column. 
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3  Main  Results 


Theorem  1  F(s)  is  robustly  -D-stable  if  and  only  if  Fe(s)  is  robustly  ZTstable. 

Proof:  Necessity  is  obvious.  To  prove  sufficiency,  we  first  note  that  interchanging  any  two  rows  (or 
columns)  does  not  affect  the  stability  of  a  polynomial  matrix  (it  only  changes  the  sign  of  the  determinant). 
By  Lemma  2 


<G>  for  all  *  =  1 


<t=>  for  all 


an(s) 


F(s)  is  robustly  D  stable 


«nl(s) 


A{8) 


an(s)  G  Kij(s) 
an(s)  G  Eij(s) 


an(s) 


ai2(s) 

:  :  ^i(s) 

«nl(s)  an2(s) 

is  robustly  D  stable. 


I  ^  i 
l  =  i 


is  robustly  D  stable. 


an(s)  G  Kij(s) 
an(s)  G  Eij(s) 
aii{s )  G  Kij(s) 
aii{s )  G  Eij(s) 


l  4  i\ 
l  =  i\ 
l  4  *2 
l  =  *2 


where  Ai(s)  is  the  corresponding  nx  (n  —  2)  polynomial  matrix.  This  last  equivalence  is  based  on 
Lemma  2  and  the  fact  that  interchanging  two  columns  does  not  change  the  stability  of  a  polynomial 
matrix.  Repeating  the  process  above,  let  Yn  denote  the  set  of  all  mappings  from  {1, . . . ,  n}  to  {1, . . . ,  n}, 
then 

F(s)  is  robustly  D  stable 

^  for  all  r 7  G  Yn,  {(««(*))  :  g  e*-(s)  i  t  v(j)  }  iS  r°bUStly  ° 

If  there  exists  an  ?y  G  Yn  such  that  rj(i i)  =  77(72)  =  k,  then  the  corresponding  matrix  F(s)  =  ( W) 
satisfies 

Oiife(s)  G  Eilk{s) 

^t2 k('^)  G  Ei2k(s) 

Applying  Lemma  2  to  column  k  of  F(s),  we  have 


F(s)  is  robustly  D  stable 
^  for  all  a  G  Sn,  {(a«(*))  :  | 

<t=>  Fe(s)  is  robustly  D  stable. 


i  ±  <r(j) 
i  =  °{j) 


is  robustly  D  stable. 


4  Interval  Model 

Interval  model,  as  a  simple  and  effective  approximation  of  uncertain  systems,  has  been  the  subject 
of  study  in  robustness  analysis  for  a  long  time.  In  a  similar  vein,  we  consider  the  Hurwitz  stability  of  the 
following  uncertain  system. 

&(s)  =  {(ct?(s))  •  cij(s )  G  Cij(s)} 

Cij  (s)  are  interval  polynomials  ' 

Definition  8  For  the  interval  polynomial  Cij(s)  =  {^D”=o  ®(*j)  G  [q,(ij),  9;(*j)]},  its 

Kharitonov  vertex  set  and  Kharitonov  edge  set  are  defined  respectively  as 

^(8)  =  {£(8), 4(8), 4(8), 4(8)} 

E'jis)  =  {Ac[.(s)  +  (1  -  A)4(s),  (r,  t)  G  {(1, 2),  (2,4),  (4, 3),  (3, 1)},  A  G  [0, 1]} 


where 

40)  =  90(*j)  +  +  92  (ij)s2  +  q3(ij)s3  +  ■■■  40)  =  qQ(ij)  +  q^ifis  +  q2(ij)s2  +  g3(ij)s 3  +  . 

40 )  =  9o  (ij)  +  g^ifis  +  g2(ij)s2  +  q3s3  +  ...  40)  =  q0(ij)  +  q^ifis  +  q.2(ij)vs2  +  q3(ij)s 3  + 
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Definition  9 


Se(s) 


u 

■es„ 


^  ( O  j  (^))nxn 


•  Cij(s ) 


e  EUs)  if  i  =  a(j)  \ 
€  Klj{s)  if  *  ^  cr{j)  J 


(4) 


Lemma  4  (Box  Theorem[16])  Suppose  A(s)  =  {<5(s,p)  =  Pi(s)Pi(s)  +  . . .  +  PTO(s)Pm(s)},  P,;(s) 
is  an  interval  polynomial,  Pj(s)  is  a  given  fixed  polynomial,  i  =  1  And  suppose  A(s)  is 

degree-invariant.  Then,  A(s)  is  Hurwitz  stable  if  and  only  if  A e(s)  is  Hurwitz  stable,  where  A e(s)  = 
UlMEtl  Fi(s)K°Pi(s)  +  ^i{s)Epi  (s)  +  YhLi+i  Fi(s)K°P%(s)}  (let  £*=r  /i  =  0,  if  r  >  t). 

By  resort  to  the  Box  Theorem,  and  following  a  similar  line  of  arguments  as  in  the  proof  of  Theorem 
1,  we  can  get  some  analogous  stability  verification  results  for  the  interval  model. 

Lemma  5  Suppose  A(s)  is  a  given  n  x  (n  —  1)  polynomial  matrix.  Then 


^-ilO)  €  f/zlO),  I 

i  =  1, . . .  ,n  [ 

f  /  cn(s)  \ 

<G>  for  all*  =  1, . . . ,  n,  <  :  A(s) 

{  \  Cnl(s)  / 


is  robustly  Hurwitz  stable. 


e  Ep(s)  l  —  i  |  S  robustly  Hurwitz  stable. 


Lemma  6  Suppose  B(s)  is  a  given  (n  —  1)  x  n  polynomial  matrix.  *  stands  for  fixed  entries  in  a 
matrix.  Then 


i(s)  *  cy(s)  * 
B(s) 

Cli(s)  *  Cij-(s) 
B(s) 


OliO)  G  f-1 i(s) 
clj(s)  €  ^1  j(s) 


is  robustly  Hurwitz  stable 


Cli(s)  x  Cij(s)  G 

(Kli(s)  x  £^0))  U  (E&s)  x  K^s)) 


is  robustly  Hurwitz  stable. 


Theorem  2  G(s)  is  robustly  Hurwitz  stable  if  and  only  if  Ge{s)  is  robustly  Hurwitz  stable. 

Remark:  Theorem  2  is  consistent  with  the  result  in  [15].  In  [15],  the  authors  obtained  their  result 
using  some  theorem  in  signal  processing.  Our  proof  is  based  on  the  properties  of  matrix  determinant, 
hence  is  more  straightforward  and  self-contained. 


5  Example 

Consider  the  uncertain  polynomial  matrix 


-40)  = 

an(s) 
012  0) 
013(5) 
021  0) 
022(5) 
023(5) 
031 (s) 
032 (s) 
033(5) 


{(a«(s)) 

=  An&ln 
=  A126I12 
=  A136I13 
=  A216I21 
=  A22W22 
=  A236I23 
=  A316I31 
=  A326I32 
=  A336I33 


3x3 


} 


0)  +  (i 
0)  +  (i 
0)  +  (i 
0)  +  (i 
0)  +  (i 
0)  +  (i 
0)  +  (i 
0)  +  (i 
0)  +  (i 


An)62n(s) 

-  Ai2)62i2(s) 
Al3)62i3(s) 
A2l)622l(s) 

-  A22) ^222 (5) 
A23)^223 (s) 

A3i)623i(s) 
A32)6232(s) 
A33)&2  33(5) 
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Then, 


Let 


Eij 

(*) 

=  Wij 

00} 

Kij 

;(S) 

=  {bl 

ij  (^)  ■  j 

(«)} 

s3 

=  {CTI, - 

,o-6} 

o-i 

1  - 

- 1;2 

~>2;3^ 

3 

cr  2 

1  - 

2;  2 

->3;3-* 

1 

03 

1  - 

->  3;  2 

— 1;3  — . 

2 

<t4 

1  - 

- 1;2 

->3;3-^ 

2 

0-5 

1  - 

->  2;  2 

— 1;3  — . 

3 

^6 

1  - 

-  3;  2 

—  2;3  — . 

1 

e  Eij(s)  if  i  =  a(j)  1 
€  Kij(s)  if  *  ±  a(j)  J 


By  Theorem  1,  A(.s)  is  robustly  D  stable  if  and  only  if  Ae(s)  is  robustly  D  stable. 


6  Conclusions 

This  paper  discussed  the  robust  D-stability  problems  for  MIMO  uncertain  systems.  The  Edge 
Theorem  and  Kharitonov  Theorem  have  been  generalized  to  multivariable  case. 
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